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TRANSFER FUNCTION DETERMINATION OF THE PRIMARY LOOP 


OF A CONCEPTUAL NUCLEAR BRAYTON SPACE POWER PLANT 
by Edward J. Petrik and Arthur W. Kieffer 
Lewis Research Center 

SUMMARY 

A phase variable technique was used to obtain transfer functions of the primary loop 
of a conceptual nuclear Brayton space powerplant operating at design conditions. 

A set of nonlinear differential equations were written describing the transient behav- 
ior of the primary loop. These nonlinear differential equations were simplified and lin- 
earized about the system design operating point. Transfer functions were obtained for 
the primary loop from the linear differential equations for small step input disturbances 
in each of the following: (1) a 1 cent step reactivity insertion, (2) a 25 percent step in 
lithium flow rate in the primary loop, (3) a 5 percent step in argon flow rate in the power 
conversion loop. 

A digital computer simulation was used to. determine the time response for the three 
input disturbances of (1) the original nonlinear differential equations, (2) the simplified 
nonlinear differential equations, (3) the linearized differential equations, (4) and the 
transfer functions. The results show that the time response obtained from the original 
equations agree with the time response obtained from the transfer functions. 


INTRODUCTION 

In the coming years, the electric power requirements for the Nation projected space 
programs will continue to increase. To meet these projected requirements, the Lewis 
Research Center has been participating in a technology program aimed at the design of a 
relatively high-powered nuclear Brayton space powerplant. 

The heat source being considered for the Brayton cycle space powerplant is a com- 
pact, fast-spectrum nuclear reactor. The design thermal output of the reactor is 2. 17 
megawatts, and the design operating lifetime for the reactor and powerplant has been set 
at 50 000 hours. 



An area of concern in the overall design of this powerplant is the design of the con- 
trols. Analysing the requirements of a closed-loop control system requires the ability 
to predict the system performance. This ability and the precision of the results dei^nd 
on how well the characteristics of each component can be expressed mathematically. 

One cannot intelligently undertake the design of any closed-loop system without a 
knowledge of the open-loop transfer function. In this report, the detailed nonlinear dif- 
ferential equations representing the primary loop of a conceptual nuclear Brayton space 
powerplant are simplified and linearized with respect to three input disturbances. The 
input disturbances considered are step changes in (1) reactivity (2) lithium flow rate in 
the primary loop (3) argon flow rate in the power conversion loop. From these simple 
linear equations, a phase variable or an indirect matrix technique (ref. 1) is used to 
transform the differential equations into phase variable form. From the phase variable 
representation, the transfer functions are directly determined. 


SYSTEM DESCRIPTION 

The equations used to simulate the primary loop of a conceptual nuclear Brayton 
space powerplant on an analog computer are presented in reference 2. In this section, 
the primary loop of a nuclear Brayton space powerplant will be briefly described and the 
equations used in the simulation will be presented. 

A simplified sketch of the primary loop of the conceptual nuclear Brayton space pow- 
erplant is shown in figure 1. The primary loop is coupled to the gas power conversion 
loop by a counterflow heat exchanger. The coolant in the primary loop is liquid lithium, 
and the working fluid in the gas power conversion loop for this study is argon. The de- 
sign point operating conditions for the primary loop of the nuclear Brayton space power- 
plant are listed in table I. 


Reactor Description 

The reactor core has 253 cylindrical fuel pins, each with a diameter of 3/4 -inch 
(1. 905 cm) and a length of 14. 8 inches (37. 6 cm). The pins are made of uranium nitride 
and are clad with a tantalum alloy with a tungsten liner between fuel and clad. The inte- 
rior of each fuel pin has a central void. The pins in the assembly are cooled by lithium 
which flows through annular flow passages formed by the outside surfaces of the fuel pins 
and the inside surfaces of the surrounding tantalum tubes. The tantalum tubes have an 
inside diameter of 0. 830 inch (2.11 cm). The reactor power is regulated by six control 
drums which, when rotated, move fuel (or poison) in or out of the core region. 
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The model used to represent the core consists of a single fuel pin with the composi- 
tion and dimensions given previously. For calculation purposes, the fuel pin model was 
divided into three axial segments of equal length. Each segment, consisting of fuel, 
cladding, the lithium coolant, is described by its proper heat transfer equation. A 
sketch of the single fuel pin model of the core is shown in figure 2. 


Core Analytical Equations 

The equations used to simulate the reactor core are described briefly in this section. 
Core kinetics . - The reactor thermal power is determined from the following 
equation: 


dt 




X.Di 


( 1 ) 


(All symbols are defined in appendix A. ) 

As indicated by equation (1), six delayed neutron groups were considered. 

The rate of change in concentration of the delayed-neutron precursors is given by 




dt 


I 


^Q-Xp. 


( 2 ) 


The decay constants and yield fractions of the precursors for the six delay groups are 

* -8 

listed in table n. The prompt neutron generation time, I , is 4. 0x10 seconds. 

Reactivity . - The reactivity feedback equations include three separate feedback ef- 
fects. They are (1) the Doppler coefficient for the reactor fuel and cladding (2) the core 
geometric expansion temperature coefficient, and (3) the density temperature coefficient 
for the lithium coolant. 

The reactivity feedbacks due to the Doppler effect, the core expansion, and the lith- 
ium density change were calculated from the following equations : 


d(6kjj) 

dT 


^ , 2,17X10-^ _ „ ^-0. 8 

dT„._ 


F, av 


d(6k ) 

- 5.69x10 


dT 


Lj « 3.V 


(3) 

(4) 
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( 5 ) 


d(6k ) 

-^= - 1.627 


dT 


L, av 


The terms Tj, Tj^ and in equations (3) and (4) represent the average 

temperatures of the fuel, cladding, and lithium in the core, respectively. These aver- 
age temperatures were determined by wei^ting the respective temperatures in each 
axial segment of the core in proportion to an assumed axial power distribution. 

The axial power distribution was obtained from a digital calculation at the design op- 
erating point. The shape of the axial power distribution at the design point resembles a 
chopped cosine. The shape of the axial distribution was assumed to remain the same at 
all thermal power levels. The power in the three axial segments of the core was distrib- 
uted such that about 39.0 percent was produced in the middle segment (segment 2) and 
about 30, 5 percent was produced in each of the other two segments (segments 1 and 3). 

Core heat -transfer . - The heat transfer equations for the core, based on a single 
fuel pin model, are presented below. 

The fuel temperature in each axial segment of the core was determined from the 
equation 


_ ^r,av^m,F ' ^K ' '^F ^ KQ 

dt (pVCp)p Arj, (pVCp)p 

and the cladding temperature was determined from the equation 

_ ^K,av^m,K '^F ~ '^K , ^L,C^C^'^L ~ 

dt (pVCp)j^ Arg. 

The lithium temperature in each axial core segment was determined from 

^'^L,C _/w\ ^'^L,C ^L,C®C - V 


( 6 ) 


(7) 


( 8 ) 


Heat Exchanger Description 

Model . - The heat exchanger model described in this section should not be consid- 
ered as a design configuration. Rather, it is only one of many possible conceptual 
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configurations capable of transferring the heat load from the primary loops to the gas- 
power conversion loop. 

A counterflow she 11-and -tube type heat exchanger was assumed for the simulation. 
The conceptual heat exchanger has 331 circular tubes; the tubes have an inside diameter 
of 3/4 -inch (1.905 cm), a wall thickness of 0.05 inch (0. 127 cm) and a length of approx- 
imately of 8 feet (2. 44 m). 

Argon gas is assumed to flow inside the tubes, and lithium flows countercurrently in 
the shell. The heat exchanger tubes are arranged in a hexagonal cross sectional array. 

The model for the counterflow heat exchanger simulation is a single tube with argon 
flowing inside the tube and lithium flowing countercurrently on the outside. The model 
was considered to be made up of three axial segments of equal length. Figure 3 is a 
sketch of the heat exchanger model. 

Heat -exchanger equations . - The heat transfer equations used in the simulation of 
the heat exchanger are described briefly below. 

The lithium temperature in each heat exchanger segment was computed from the 
equation 


dT 


L, HX _ ^'^L, HX , ^L, HX^L, HX 

dt p^Al AXjjx PlJlfp,L 


HX ■ HX^ 


(9) 


And the wall temperature for each segment of the heat exchanger was computed from 
the equation 


HX _ ^L, HX^L, HX 
dt ■ (PVCp)^ 


h s. 


^■^L,HX ■ '^W,HX^ + ^^pV " ^W,HX^ 


( 10 ) 


The temperature of the argon gas in each axial segment was computed from the 
equation 


dt p 


g 


^g ^'^g,HX , Vg,HX (^ 

^HX Vg^p, g 


W,HX 


- T ) 

cr' 


( 11 ) 


Heat -Transfer Coefficients 

The liquid metal heat transfer coefficient for the reactor core was calculated from 
the equation 
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NUl c = 5. 591 + 4. 595xlO“^^eL q)®’ (12) 

And the liquid metal heat transfer coefficient for the heat exchanger was calculated 
from the equation 



The film coefficient on the gas side of the heat exchanger was calculated from the 
Dittus-Boelter expression: 

NUg = 0. 023^Reg^°* ^ (14) 

Transport Lag Representation 

In the primary loop of the nuclear Brayton space powerplant, the transport lag (or 
time delay) is defined as the average time required for the coolant to traverse the loop. 
However, it is assumed that the heat exchanger is close coupled to the reactor such that 
nearly all of the time delay occurs in the heat exchanger. Because of this, only the 
transport lag in the heat exchanger was considered. Thus, the fluid temperatures at the 
inlet of each of the three heat exchanger segments were delayed by one -third of the delay 
time before entering the next heat exchanger segment. 

At the design lithium flow rate, the time delay is about 10 seconds. Therefore, a 
change in the lithium temperature at the reactor outlet goes unnoticed at the reactor in- 
let for about 10 seconds. For lithium flo.. rates of 50 percent and 10 percent of the de- 
sign value, the time delays would be about 20 seconds and 100 seconds, respectively. 

It should be recognized that these time delays are primarily a characteristic of the 
heat exchanger design used in the simulation. These time delays may change with the 
heat exchanger design configuration. 

ANALYSIS 

Analysis of a dynamic system requires the ability to predict the systems perform- 
ance. The ability and precision of the results depend on how well the characteristics of 
each component can be expressed mathematically . In the previous section, the equations 
(ref. 2) used to simulate the nuclear Brayton space powerplant on an analog computer 
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were presented. In this section the nonlinear differential equations will be simplified 
and then linearized about the design operating point of the nuclear Brayton space power- 
plant (NBSP). Transfer functions will be derived from the linear differential equations 
for small step input disturbances. 


Simplifying the Brayton System Equations 

In general, three simplifying assumptions will be used to reduce the number of non- 
linear differential equations which describe the NBSP. 

The first assumption considered is related to equation (1). As indicated, six delayed 
neutron groups were considered in the original simulation. It will be assumed that there 
is only one group of delayed neutrons. And the decay constant is equal to a properly 
weighted average of the six actual groups (ref. 3). 

The decay constant is computed as follows 

X= = 7.85x 10"2 sec"^ (15) 

6 

^i 



The next simplification considers one axial segment (lump) for both the reactor and 
heat exchanger. Originally, each model was divided into three axial segments or lumps. 

Because the reactor heat transfer equations will now consist of one axial segment, 
the calculated fuel and cladding temperature will represent the average fuel and cladding 
temperature. Thus, the terms T^ and T^- in the reactivity feedback (eqs. (3) 
and (4)) will be Tj, and Tj^. In addition, the average temperature of the lithium in the 
core for the reactivity feedback will be determined by averaging the reactor inlet and re- 
actor exit lithium temperature. This is expressed by equation (16). 

T ^ '^ L>out '' jL,in (16) 

JLi f 3.V n 


Finally, it is assumed that the reactor cladding temperature is equal to the lithium 
temperature in the core. And the wall temperature of the heat exchanger is equal to the 
lithium temperature in the heat exchanger. This assumption is valid for step disturb- 
ances in reactivity and argon flow rate since the resistance to heat transfer between the 
cladding and liquid or liquid to heat exchanger wall is negligible . That is, the liquid 
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metal heat transfer coefficient is extremely large. Therefore, an overall heat transfer 
coefficient is used in the simplified nonlinear differential equations. 

For a step disturbance of lithium flow rate the preceding assumption that the reactor 
cladding temperature is equal to the liquid lithium temperature in the core is valid. How- 
ever, the wall temperature associated with the lithium temperature at the heat exchanger 
exit is not negligible. Thus, the heat exchanger wall temperature will be included in the 
simplified nonlinear differential equations. 

Applying the above simplifications to the original nonlinear differential equations re- 
duces the number of equations from 25 to 5 for step disturbances of reactivity and argon 
flow rates . And reduces the number of equations from 25 to 6 for a step disturbance of 
lithium flow rate. Appendix B presents the final simplified nonlinear equations. 


Linear Brayton System Equations 


The set of simultaneous nonlinear differential equations as derived are linearized in 
this section. These equations are expressed as five first-order equations of the form 


^ = f (X, u) (17) 

dt 


where 


u 


is the input. 


The shorthand notation for the column vector is a 


as example 



X = 


= col(Xj,X2, 


V 



( 18 ) 


To find a linear approximation of the nonlinear equations, the equations must be linear- 
ized about some operating point. The solutions are then valid only in a neighborhood of 
the point. 

^ ^ A 

The equations are linearized by expanding f(X, u) in a Taylor series (ref. 4) about 
the operating point X° and u°, as follows 


f (X, u) = f(i^, u°) + 

ax 


(X - X°) + 

X=^ 3u 


u=u 


(u-u°) + higher -order terms (19) 

x=5^ 


u=u 
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Let X = X - and u^ = u - u be the changes in X and u from the operating point 

X° and u°. Since X° and u° are assumed to be equilibrium conditions, 

dX/dt = dX*/dt and f(X°,u°) is defined to be 0. Thus, the first-order approximation of 

^ = f (X, u) (20) 

dt 


can be written as 


dX* _ 8f (X, u) 

dt ax 




3f(X,u) 

3u 



u=u 


u=u 


o 


The linearized set of differential equations now take the form 


dX* 

dt 


= + Bu* 


( 21 ) 


( 22 ) 


where A and B are matrices, the elements of which are determined from 




3f.(X,i 

A A 

x=x' 

3x. . 

A Ar 

u=u 


and 


3f.(X,u) 

b.. = — 

13 


3u. 


|x=x° 

u=u 


(23) 


(24) 


The linearization procedure as described above assumes that the vector function is 
continuously differentiable in all its arguments . Two variables, the transport delay and 
the heat transfer coefficients on the gas side of the heat exchanger, are not continuously 
differentiable . 

The transport delay will be approximated by using a second-order Fade approxima- 
tion given in reference 5. A higher order approximation could be used except that the 
number of differential equations needed for the approximation is equal to the order of the 
approximation. Thus, the transport delay y = x where x is delayed by t seconds is 
given as follows. 
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(25) 


du 

dt 


V 


12 


X 


T 


dv 

dt 


12 6 72 

— u — v + — X 


y = u + X 


(26) 


(27) 


The heat transfer coefficient is a discrete function of the gas temperature given in 
tabular form. A least-square quadric curve fit was used to make the heat transfer coef- 
ficient a continuous function of gas temperature without increasing the number of differ- 
ential equations. 

Appendix C presents the detailed linearization as described previously for reactivity 
and argon flow rate disturbances. The linearization procedure for lithium flow rate dis- 
turbances is the same as appendix C except the nonlinear differential equations include 
the heat exchanger wall temperature. 


Transfer Functions 

If a system is linear and time invariant, the transfer function conveniently relates 
the response of the system to a forcing function with the aid of appropriate Laplace trans- 
forms. If R(S) represents the Laplace transform of the input disturbance, C(S) repre- 
sents the transform of the response, and G(S) represents the transfer function, then 

C(S) = R(S)G(S) (28) 


or 


input 

R(S) 


G(S) 


output 

C(S) 


In the preceding equation or diagram it is assumed that the system is initially at rest, 
that is, all initial conditions are equal to zero. 

Since the equations representing the nuclear Brayton space powerplant derived in the 
previous section are linear differential equations with constant coefficients, the basic 
problem associated with determining the transfer functions is solving the equations only 
in terms of the output and input variable. This becomes tedious if the order is greater 
than two. 
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Block diagram manipulation can also be used to find the transfer functions. Unfor- 
tunately, this also becomes laborious if there are many internal feedback loops such as 
the reactivity feedback in the NBSP. 

The method used in this report is referred to as the indirect matrix or phase var- 
iable approach. This method is based on the realization that, if the system is trans- 
formed into phase variables, the transfer functions may be determined by inspection. In 
addition to determining an overall transfer function, any number of internal transfer func- 
tion may be determined at the same time. Reference 1 formulates the indirect matrix 
method and describes a digital program implementing the method. 

The indirect matrix or phase variable method was applied to the linear differential 
equation to obtain the transfer functions. The list of transfer functions in both polynomial 
and factored forms of the Laplace transform variable s for the three input disturbance 
follow: 

Input disturbance in reactivity . - 
Denominator: 

D(s) = s'^ + 33. 3 s® + 224. 5 s^ + 224. 7 s^ + 102. 5 s^ + 27. 0 s^ + 0. 494 s + 0. 01493 


or 

SpS2 = - 8.653X10"^ ± 2.277xl0"^ j 
S3, s^ = - 0. 2288 ± 0. 4005 j 
S5 = - 0.6279 
Sg = - 7. 688 
s^ = - 24.51 


Numerator: 

(1) Concentration of delayed-neutron precursors 

C(s) 

Pjj^(s) D(s) 

N^(s) = [1. 341 s® + 44. 67 s^ + 300. 8 s^ + 294. 0 s^ + 130. 4 s^ + 34. 65 s + 0. 3239]xl0^ 
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or 


Sj, S 2 = - 0.2151 ± 0.4016 j 
S 3 = - 0.6344 
84 = - 0.9693X10"^ 

S 5 = - 7. 718 
Sg = - 24. 51 

(2) Reactor fuel temperature 

Tji(s) N2 (s) 

Pjj^(s) D(s) 

N 2 (s) = 0. 8032 s® + 26. 13 s® + 162. 5 + 115. 7 + 46. 15 + 4. 811 s + 0. 1409 


or 


Sj,S 2 = - 0.3096 ± 0.3704 j 

5 3 = - 0.7854x10"^ 

5 4 = - 0.5383x10“^ 

Sg = - 7.266 

Sg = - 24.51 

(3) Reactor exit lithium temperature 

Pin(s) D(s) 

NgCs) = 2. 963 s^ + 75. 13 s*^ + 62. 43 s^ + 20. 63 + 2. 838 s + 0. 1256 
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or 


SpS 2 = - 0.2879 ± 0. 1647 j 
Sg = - 0.8303x10"^ 

S 4 = - 0. 1893 
Sg = - 24.51 

(4) Heat exchanger exit lithium temperature 

p.j^(s) D(s) 

N^Cs) = 0. 4969 + 11. 94 - 5. 972 + 0. 7833 s + 0. 1041 

or 


Sp S 2 = 0. 2857 ± 0. 1650 j 

5 3 = - 0. 7854x10'^ 

5 4 = - 24.51 

(5) Argon gas temperature 

Tg(s) _ Ng(s) 

Pin^s) D(s) 

Ng(s) = 57. 21 + 48. 53 + 16. 17 + 2. 232 s + 0. 9694x10“^ 


or 


Sj, S 2 = - 0.2857 ± 0. 1650 j 
Sg = - 0. 1982 
S 4 = - 0.7854X10'^ 
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Input disturbance in argon flow rate . - 
Denominator: 

D(s) = + 34. 47 s® + 234. 7 + 235. 2 + 107. 2 + 28. 17 + 0. 5262 s + 0. 1554x10"^ 

or 

Sj, S2 = - 0. 8875x10"^ ± 0. 2267x10"^ j 
Sg, = - 0.2292 ± 0. 3991 j 
Sg = - 0. 6270 
Sg = - 7.687 
= - 25. 68 

Numerator: 

(1) Concentration of delayed -neutron precursors 

C(s) _ 

Wg(s) D(s) 

Nj(s) = [12.78 s^ + 499.2 s^ + 4573.0 + 4793.0 s^ + 1774.0 s + 248,9]xl0“'^ 

or 

SpSg = - 0.2857 ± 0. 1650 j 
Sg = - 0.5697 
s^ = - 12.23 
Sg = - 25. 68 

(2) Reactor fuel temperature 

Tp(s) NgCs) 

Wg(s) " d”(s) 
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NgCs) = 7. 655 - 217. 4 s'* - 0. 1086x10^ - 5973. 0 - 1063. 0 s + 20. 07 

or 

Sp S2 = - 0.2857 ± 0. 1650 j 
S3 = 0.1718X10"* 
s^ = - 25. 68 
Sg = 54.63 

(3) Reactor exit lithium temperature 

Wg(s) ~D(s) 

NgCs) = - 598.4 s^ - 0. 1619x10^ s^ - 0.2148x10^ s^ - 8729.0 s^ - 1312.0 s + 6.348 
or 


Sp Sg = - 0.2857 ± 0. 1650 j 

53 = 0.4690x10“^ 

54 = - 0.8093 
Sg = - 25. 68 

(4) Heat exchanger exit lithium temperature 

'*'l, 

Wg(s) D(s) 

N^(s) = -[171. 4 s® + 5891. 0 s® + 0. 3928x10® s'* + 0. 3524x10® s® + 0. 1203x 10® s^ + 


1574.0 s + 8. 198] 
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or 


Sp Sg = - 0.2857 ± 0. 1650 j 
Sg = - 0.5430x10"^ 

S 4 = - 0.4134 
Sg = - 7.602 
- --25.68 

(5) Argon gas temperature 

Tg(s) _ Ng(s) 

Wg(s)~D(s) 

Ng(s) = -[16. 82 s® + 147.9 s^ + 163.4 s^ + 86. 03 s^ + 25. 23 + 1. 37 s + 0. 5216x10"^] 10^ 

or 

SpSg = - 0.2569 ± 0.4370 j 

5 3 = - 0.4111X10'^ 

5 4 = - 0.6308X10'^ 

Sg = - 0.6123 

Sg = - 7.597 

Input disturbance in lithium flow rate . - 
Denominator: 

D(s) a s® + 56. 17 s"^ + 1001.0 s® + 5759.0 s^ + 6317.0 s“^ + 3152.0 

+ 804.0 + 12.27 + 0.3609 
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or 


Sj, S 2 = - 0, 7139x10"^ ± 0, 2063x10"^ J 
S 3 , = - 0.2940 ± 0.4125 j 

Sg = - 0.6473 
Sg = - 8.465 

Srj = - 22.28 

Sg = - 24.17 

Numerator: 

(1) Concentration of delayed -neutron precursors 

C(s) 

Wl(s) D(s) 

Nj(s) = [4. 235 s® + 203. 0 s^ + 2574. 0 + 3460. + 1324 s^ + 43. 79 s - 0. 0289]xl0® 

or 

Sj,S 2 = - 0.7145 ± 0. 1242 j 

5 3 = 0.6471x10"^ 

5 4 = - 0.3718x10'^ 

Sg = - 22.31 

Sg = - 24. 17 

(2) Reactor fuel temperature 

Tp(s) _ N 2 (s) 

Wl (sj ~ dST 
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N 2 (s) = -[0. 2091x10"^ + 3. 474 s® + 160. 3 + 1931 + 1342 

+ 339. 0 + 4. 525 s + 0. 151]10^ 

or 

s j, S 2 = - 0. 6064x10“^ ± 0. 2079x10“^ j 
83,84 = - 0.3553 ± 0.2261 j 

Sg = - 22.26 

8g = - 24. 17 
8^ = - 1614.0 

(3) Reactor exit lithium temperature 

^L, ^3^®^ 

Wl(s) D(^ 

N 3 ( 8 ) = -[4. 23 + 203.2 8 ® + 2590. 8 ^ + 3745. 8 ^ + 1843. 

+ 352.7 8 ^ + 4.911 8 + 0. 1500]xl0^ 
or 

8 J, 82 = ~ 0. 6224x10”^ ± 0. 2043x10"^ j 
83,84 = - 0.3571 ± 0.2062 j 
85 = - 0.8486 

Sg = - 22.26 

8^ = - 24. 17 

(4) Heat exchanger exit lithium temperature 

'^L, _ ^4^®^ 

Wl(s) D( 8 ) 
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N^Cs) = [0.1152 8”^+ 5.811 s® + 90.08 + 437.8 s'^ + 715.9 

+ 416. 6 + 6. 664 s + 0. 2020 ]k 10^ 

or 

®1’ ®2 " ■ '^'i'SSxlO"^ ± 0.2092 x 10"^ j 
83,8^ = - 1. 188 ± 0.4929 j 
Sg = - 4. 555 
Sg = - 19.33 
8,^ = - 24.17 

(5) Argon gas temperature 

Tg(s) _N5 (s) 

^l(s) D(^ 

Ng(8) = [3. 103 8"^ + 17. 85 8® - 1234. - 2126. - 1055. 

- 177.4 - 2 . 330 s - 0. 7302xl0"^]10^ 

or 

Sj, 83 = - 0. 5685X10“^ ± 0. 2027x10"^ j 
S3, 8^ = - 0. 3465 ± 0. 1082 j 
Sg = - 0.9938 

8g = 18.21 

8 y = - 22.27 
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Of special significance is the order of the characteristic equations for the various 
disturbances. For a lithium flow rate disturbance, the characteristic equation is eighth- 
order whereas the characteristic equation is seventh -order for reactivity and argon flow 
rate disturbance. 


RESULTS 

In this section, the time response of the reactor exit lithium temperature to a step 
disturbance in (1) reactivity, (2) argon flow rate, and (3) lithium flow rate will be de- 
scribed. The time response was determined by digital computer using the continuous 
system modeling program. The step disturbances were imposed with the nuclear Brayton 
space powerplant operating at design (table I). 

As stated in the ANALYSIS section, the original nonlinear differential equations 
(ref. 1) were simplified and linearized before obtaining the transfer functions. To show 
the effects of the assumptions, the time response of the reactor exit lithium temperature 
will be presented at each step of the simplification. The assumptions used are described 
and referred to as case A to D in table III. 


Step Change in Reactivity 

With the nuclear Brayton space powerplant at the steady -state design point, a pos- 
itive step change in reactivity of one cent was made. As a result, the thermal power out- 
put of the reactor changes and the lithium temperature at the reactor exit changes with 
time as shown in figure 4. 

In what follows, the difference in time response of the reactor exit lithium tempera- 
ture will be described by pointing out the effects of the simplifying assumptions. 

Figure 4(a) is a plot of the time response for the original simulation. A peak tem- 
perature change of about 9. 5 °R (5. 3 K) is reached in approximately 125 seconds. Figure 
4(b) is the time response of the reactor exit lithium temperature for case A. In this fig- 
ure, a peak temperature change of approximately 11° R (6. 1 K) is reached in about 150 
seconds. The difference in the time response can be minimized by using some other ap- 
proximate formulation of the six -group delayed neutron. However, for a closed-loop 
control this difference may be considered negligible. 

Cases B and C are now considered by comparing figures 4(c) and (d) with 4(a). The 
difference in the time at which the peak temperature occurs as compared with figure 4(a) 
is less than 10 seconds. In addition, the difference in the peak value is less than 2° R 
(1.1 K). The difference in the peak value and the time at which it occurred can be con- 
sidered negligible. 
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Figures 4(e) and (f) are plots of the time response of the reactor exit lithium tem- 
perature for the linear differential equations and its corresponding transfer function. 
Comparing these figures with figure 4(a) shows that the- transfer function approximates 
the reactor exit lithium temperature. 


Step Change in Argon Flow Rate 

With the system operating at design conditions, a 5 -percent step increase is made 
in the argon flow rate. This increase in argon flow rate directly affects the system in 
two ways. It reduces the dwell time of the gas in the heat exchanger and it also improves 
the gas heat transfer coefficient. 

Figure 5(a) is a plot of the lithium temperature at the reactor exit for the original 
simulation. The reactor exit temperature decreases during the initial 50 seconds be- 
cause more heat is transferred to the gas in the heat exchanger. 

Figures 5(b) to (d) are plots of the time response of the reactor exit lithium temper- 
ature for case A to C. Comparing these figures with figure 5(a) also shows that the sim- 
plifying assumptions are valid as they were for the reactivity disturbance described pre- 
viously. In addition, the transfer function approximates the reactor exit lithium temper- 
ature for a step disturbance in argon flow rate since figures 5(e) and (f) compares with 
figure 5(a). 


Step Change in Lithium Flow Rate 

With the NBSP at steady-state design conditions, a step increase of 25 percent was 
made in the lithium flow rate. A large step in lithium flow rate was chosen because the 
liquid lithium flow rate effects on the NBSP are rather insensitive. This behavior is ex- 
plained in detail in reference 2. 

The time response of the reactor exit lithium temperature to the step in lithium flow 
rate for the original simulation is shown in figure 6(a). 

The immediate effect of this disturbance is an increase in the lithium velocity 
through the reactor and heat exchanger. This velocity increase has little effect on the 
lithium heat transfer coefficient. However, it does reduce the lithium dwell time in the 
reactor and heat exchanger. As a result, the lithium loses less heat in the heat ex- 
changer and also gains less heat in the reactor. Therefore, the lithium temperature at 
the reactor exit immediately decreases approximately 14° R (7. 8 K) during the initial 3 
seconds as shown in figure 6(a). 
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The time response of the reactor exit lithium temperature for cases A, B, and D 
are shown in figures 6(b) to (d). The initial peak occurs during the first 3 seconds (the 
same as in fig. 6(a)). And the peak values for figures 6(a) to (d) are within 2° R (1. 1 K). 

Figure 6(e) and (f) are plots of the time response of the reactor exit lithium temper- 
ature for the linear differential equations and its corresponding transfer function. Com- 
paring these figures with figure 6(a) again shows that the transfer function approximates 
the reactor exit lithium temperature of the original simulation. 

In the foregoing discussion, only the time response of the reactor exit lithium tem- 
perature to an input disturbance was described. This showed the step by step procedure 
that was used to obtain the different transfer functions which relate the reactor exit lith- 
ium temperature to the different input disturbances. In table IV, we have summarized 
the time response of the major variables comparing the original simulation and the trans- 
fer function approximation. 

Of special significance in table IV is the time of the peak overshoot. With the NBSP 
initially at its steady-state design point, the time of peak overshoot for a reactivity step 
insertion of 1 cent is about 125 seconds. For a 5-percent step in argon flow rate, the 
time of peak overshoot is about 70 seconds. In contrast, a 25 -percent step disturbance 
in lithium flow rate, the time of peak overshoot is less than 12 seconds. Hence, the re- 
sponse of the system to an input disturbance in flow rate is extremely fast. 


SUMMARY OF RESULTS 

The transfer functions representing the primary loop of a conceptual nuclear Brayton 
space powerplant operating at design power were determined for three input disturbances. 
The input disturbances included (1) a 1 cent step of reactivity, (2) a 5 -percent step of 
argon flow rate in power conversion loop, and (3) a 25-percent step of lithium flow rate 
in primary loop. 

The primary loop was described by a set of nonlinear differential equations. These 
equations were simplified and linearized about the design power operating point. The 
transfer functions were obtained from the linear differential equations by usir^ an indi- 
rect matrix or phase variable method. 
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A digital computer simulation was used to determine the time responses of the orig- 
inal nonlinear differential equations and of the transfer functions . The results indicate 
that the transfer functions do represent the primary loop .of a conceptual nuclear Brayton 
space powerplant. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 29, 1970, 

120-27. 
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APPENDIX A 



SYMBOLS 


A 

flow area, ft^; 

Nu 

Nusselt number 

A 

coefficient matrix 

Pe 

Peclet number 

An, 

mean heat conduction area, 

Pr 

Prandtl number 


fr; m"' 

Q 

reactor power, Btu/sec; kW 


matrix elements 

Re 

Reynolds number 

B 

column matrix 

Ar 

incremental distance in a direc- 


elements of column matrix 


tion radial to flow, ft; m 

C 

concentrations of delayed neutron 

S 

2 2 

heat transfer area, ft ; m 


precursors in the i*"^ group 

s 

Laplace transform variable 

S 

specific heat at constant pressure, 

T 

temperature, °R; K 


Btu/(lbm)(°R); (kW)(sec)/(Rg)(K) 
hydraulic diameter, ft; m 

T 

average section temperature, 
°R; K 


contribution to power from i*’” 

AT 

temperature increment, °R; K 


group of delayed neutron pre- 

t or T 

time, sec 


cursors 

U 

overall heat transfer coefficient. 

d 

outside diameter of tubes in heat 


Btu/(sec)(ft^)(°R); kW/(m^)(K) 


exchanger, ft; m 

u 

independent variable 

^C 

proportionality constant, Btu(cm )/ 
neutrons sec; (kW)(cm^)/ 



u 

column vector 


neutrons 


elements of column vector 

f() 

matrix function of 

V 

volume, ft^; m^ 

h 

heat transfer coefficient, Btu/ 

V 

independent variable 


(sec)(ft2)(°R); kW/(m^)(K) 

W 

flow rate, Ibm/sec; kg/sec 

K 

fraction of total power generated 

X 

column vector 

k 

in fuel section 

thermal conductivity, Btu/(ft) 

AX 

incremental distance in a direc- 
tion axial to flow, ft; m 

L 

Z* 

(sec)(°R); kW/(m)(K) 
length, parallel to flow, ft; m 
prompt neutron generation time. 

X 

^i 

independent variable 
elements of column vector 


sec 

y 

dependent variable 
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/3 

total delayed neutron precursor 

F 

fuel 


yield 

g 

argon gas 

^i 

delayed neutron precursor yield 
fraction in i^ group 

HX 

heat exchanger 


in 

inlet or input 

X 

average delayed neutron precur- 

K 

fuel cladding 


sor decay constant 

X. 

1 

delayed neutron precursor decay 

L 

liquid lithium 

constant in i^^ group 

W 

tube wall in heat exchanger 

6k 

reactivity 

P 

pertaining to the reactivity effect 

P 

density, Ibm/ft^; kg/m^ 


due to lithium density change 

Pi 

•fVi 

constant 

T 

delayed t seconds 

Pin 

reactivity insertion 

Superscripts: 

Subscripts: 

0 

steady -state operating point 

av 

average 

* 

change in variable 


C core 

D pertaining to the Doppler re- 

activity effect 

E pertaining to the reactivity effect 

due to core expansion 
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APPENDIX B 


SIMPLIFIED NONLINEAR DIFFERENTIAL EQUATIONS DESCRIBING 
THE NUCLEAR BRAYTON SPACE POWERPLANT 


Equations for Reactivity and Argon Flow Rate Disturbances 


Concentration of delayed -neutron precursors: 


dt 


i-Q -XD 
I* 


where 


Q = 


Z*XD 
P - 6k 


Reactor fuel temperature: 


dTj, 

dt 


Sc 


-Tj,) + 


Q 

(pVCp)p 


Reactor exit lithium temperature: 


dT 


(T, 

Wl 


dt \pvK l,hx-l,c^--^^Ul,c(Tf-Tl,c) 


- Tj.. r) + 

Heat exchanger exit lithium temperature: 


'p^L 


‘*'^L,HX_/w\ ^'^L,HX ®L,HX , . 

UL,HX^^g,in" L,HX^ 


dt 


3A/jjx ^^hx ^^^p^L 
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and 


AT. 


AX 


XjjHX JC (T' T' ^ 

- - ^K,L^^L,C,t " ^L,HX^ 


HX 


Argon gas temperature: 





W„ AT„ 
g g>HX ^ 


AX. 


HX 






and 


AT 


g, HX ~ .p- /rp 

^K,g^^g,in 


AX 


HX 



Equations for Lithium Flow Rate Disturbance 


Concentration of delayed-neutron precursors: 


® -i. 

dt ~ I* 


Q 


- XD 


where 


Q = 


Z*XD 
l3 - 6k 


Reactor fuel temperature: 


dTp 

dt 




-Tp) + 


Q 

(PVCp)p 


Reactor exit lithium temperature: 


dT 


L,C _ 


dt 


(s). 


(T 


L,HX " ^L,C 


) + 


(pVCp)^ 


(T, 


'^L,C^ 
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Heat exchanger exit lithium temperature: 


and 


dT 


L,HX 


dt 


w\ ^^L,HX . ®L,HX ^ frr 
+ 7 - "l,hx^^w 


" HX^ 


AT 


L , HX /rp rp > 

T^;: ^k,l'^^l,c,t ^l,hx^ 

AArx 


Argon gas temperature: 


dT„ 
g _ 

1 


^'^g,HX , 

(s \ 

U tT, ^ 

- T ) 

dt 

^g 


^^HX 

\VC / 

\ p/g 

V L,C 



and 


AT 


AX 


g, HX ^ j, /„ _ T 1 

^K,g^^g,in ^g^ 


HX 


Heat exchanger wall temperature; 


'■Tw_®L,HX ^ ) + ®g.HX ^ , , 

« fcVCJ,„ 


p'W 


28 



APPENDIX C 


LINEARIZATION OF THE SIMPLIFIED NONLINEAR DIFFERENTIAL EQUATIONS 
FOR REACTIVITY AND ARGON FLOW RATE DISTURBANCES 


The following equations were presented in appendix B, except that the constant coef- 
ficients are now expressed numerically and the variables are expressed in terms of the 
column vector elements. 

Concentration of delayed-neutron precursors: 

Xj = Xj3(j8 - 6k)"^ Xj - Xxj 

Reactor fuel temperature: 

X2 = 2. 398 Uj^ Q(xg - Xg) + 9. 88 XZ*(/3 - 6k) Xj 
Reactor exit lithium temperature: 

Xg = 42. 5 W(X^ - Xg) + 10. 37 Uj^ ^^(Xg - Xg) 

Heat exchanger lithium temperature: 


^ ^K, L^HX^^S’ " ^ 4 ^ + 2.645 jjx^'^g, : 


Second-order temperature time delay: 


where 


^3,t = ^3 + ^5 
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Argon gas temperature : 


^7 = 


2.88 


l-258%,gWg(Tg,i„-x^)+1.99Ug(x3 


The preceding seven equations can be written as follows, 


X = f(X, u) 



Define 


Then 


Pj = 2. 17x10"^ 

P 2 = - 3.233x10"^ 

Pg = - 7.317x10"® 
p^ = 1.96x10"^ 

Pjj^ = reactivity insertion 


86k _ ^1 
3 X 2 ''2 


38k 

3Xo 


= 0.2 


- 0 . 8 

P2X3 + 


2 


36k _ ^3 
3 x 4 2 

From equation (Bl), the elements of the first row of the coefficient matrix is determined. 
Note that all partial derivatives are evaluated at the design operating point (i.e. , X = X° 
and u = u°) 


3xj /3 - 6k /3 
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aj2 - 


^^1 ^ / 96k\ _ ^^1 ^1 


9Xr 


= -X|3Xj(j3 - 6k)' 


9Xr 




dXt Xx 


o 


^13 - 


1 '^^1 96k 


Xx^ 


9x, 


/3 9Xo /3 


0.2 P2(x°)“°-^ + ^ 


aj4 - 


^^1 _ ^^1 96 k ^ ^ 

9x. ^ 9x. /3 2 


^15 " ^16 ■ Hi 


9Xf 


= 0 


The remaining elements of the coefficient matrix are determined in the same manner 
as described 


9x„ 3. , * 

a„. =_i= 9.88 2^ 
9x. ^ 


H 2 " 


dx^ 

9Xr 


X?V 

-2.398Ux p + 9. 88 

3X2 


9Xr 


a„ = 2.398 U, „ + 

8x3 ^ 8x3 


a. 
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0 x 2 9.88XZ x° ggj^ 


9X/ 




9x, 


^25 “ ^26 ~ Hl~^ 


agi - 0 


as2 - 10.37 Ul c 


agg = -42.5 W - 10. 37 Ul q 


ag4= 42. 5 W 
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^35 “ ^36 “ ^37 " ° 


^41 ~H2~^ 


^43 " %, L^HX 


9x4 

a^4 = — = - 11. 8 lWjjx - 2. 645 U^ 


^45 " ^ ~ ® ^K, L^HX 


^46 - a4-j - 0 


^51 " ^52 “ ^54 " ^57 


3X| 


^53 “ 


5 _ -12 

8X3 T 


^55 " 


-12 


^56 “ 


-6 


^61 ~ ^62 “ ^64 " ^67 “ ® 


^63 " 


72 


^65 " 


■12 


^66 “ 


-6 


^71 “ ^72 “ ^74 ~ ^75 “ ^76 " ® 


,HX 
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^77 “ 


3X7 "^7 


dXrj 2 . 88 


3 x 3 2.88 S' 


^ ^ 3u„ 

- 1. 258 W„ - 1. 99 U + 1. 99 x2 - x2 — i 
K,g g g 3 7 g 


1 . 258 gWg(Tg_ in - X?) + 1. 99 Ug(a° - x?) 

Vg ^^-1. 25x10"^ x^ + 5. 75x10"® x.^ + 1. 45x10"^^ 


2.88 


h„ 3h„ U? 3h 

g^ g i. = _i i 

^^7 + 0. 46^^ ^^7 h| ^^7 

-^ = 4. 08 W®- ®f-2. 5X10"^ x2 + 5. 75x10"®) 
dXn S \ / / 


The elements of the input matrix are determined in the same manner as the coefficient 
matrix. 

Reactivity input disturbance: 


3xi Xx 


bi = 


1 '^^1 36k ''^^1 


Xx? 


9 Pin ^ 9 Pin ^ 


•^2 “ 


3X2 ^*®®^^*^136k 9-88XZ*xJ 

3pin ^2 3p.^ ^2 


ba - b4 - b5 - bg - b7 - 0 


Argon flow rate disturbance : 


bf - b2 - bs - b5 - bg - 0 
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dx. 

b4 = -A = 11.8 

^ 9W 


f o o o\ 

Xg + X3 - X^j 


9W 9W, 


^ = -0.5443 W„+ 0.2817 

K.J g 


Ul,hx = (1- + 0.46V* 

"l,hx 


h„ = 8.33x10"^ W°-® 


g,in 


g 


9K 


J^= - 0.5443 

awg 

= u; 


^“l,HX _„2 1.13 in 


9W„ .2 9W„ 

g g, in g 


= 6. 66x10"^ W"®‘ ^ 
g 


9W 


g 


^7 " 


9Xr^ Xi^ 


9W„ 2.88 


1 . 258{Tg_ - x“ L ^ + 1. 99(x» - x») 


9 U„ 

g 

3W„ 


Kt;. „ = 0.219 W„+ 0.0366 
g g 


D^=(h-1.0.46)-l 


h„ = 4.08 W 


. 0 . 3 (.,. 


25X10"^ + 5. 75x10"® x,y + 1. 45x10"^ 


9K 


^ = 0.219 

9 W„ 
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au 

aw 


g A 




3.26W‘°'^ (-1.25x10"® x2-x2+ 5,75'® x2 + 1.45X10'®) 

aw ^ ^ 

o 
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TABLE I. - DESIGN POINT OPERATING CONDITIONS 


Reactor thermal power, Q, MW^, 

2. 17 

Lithium temperature at reactor inlet (or heat 

2100 (1167) 

exchanger exit), °R (K) 


Lithium temperature at reactor exit (or heat 

2200 (1222) 

exchanger inlet), °R (K) 


Lithium flow rate, Wj^, Ibm/sec (kg/sec) 

20.7 (9.39) 

Argon temperature at heat exchanger inlet, °R (K) 

1560 (867) 

Argon temperature at heat exchanger exit, °R (K) 

2060 (1144) 

Argon flow rate, Wg, Ibm/sec 0^/sec) 

32.8 (14.88) 
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TABLE n. - DECAY CONSTANTS AND YIELD FRACTIONS 


OF DELAYED-NEUTRON PRECURSORS - 


FAST FISSION 


Group 

Decay constant, 
Xi, 

sec'^ 

Yield fraction, 

1 

3. 8800 

0.0001719 

2 

1.4000 

. 0008430 

3 

. 3110 

. 0026900 

4 

.1160 

. 0012380 

5 

.0317 

. 0014020 

6 

.0127 

. 0002520 


TABLE m. - NUCLEAR BRAYTON SPACE POWERPLANT DESCRIPTION AND SIMPLIFYING 
ASSUMPTIONS OF THE NONLINEAR DIFFERENTIAL EQUATIONS 


Case 

Reactor kinetics 

Reactor heat transfer 

Heat exchanger heat transfer 

Original 

simulation 

Six neutron groups 

Three axial segments in- 
cluding cladd, fuel, and 
lithium temperatures 

Three axial segments in- 
cluding argon gas, wall, and 
lithium temperatures 

Case A 

- One delayed neutron 
group 

Same as original simulation 

Same as original simulation 

Case B 

One delayed neutron 
group 

One axial segment including 
cladd, fuel, and lithium 
temperatures 

One axial segment including 
argon gas, wall, and lithium 
temperatures 

Case C 

One delayed neutron 
group 

One axial segment including 
fuel and lithium tempera- 
tures, i. e., temperature of 
lithium to clad assumed zero 

One axial segment including 
argon gas, and lithium tem- 
perature, i. e., temperature 
of lithium to wall assumed zero 

Case D 

One delayed neutron 
group 

Same as case C 

Same as case B 
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TABLE IV. - SUMMARY OF TIME RESPONSE FOR VARIOUS STEP DISTURBANCES 



Original simulation 


Peak over- 
shoot 





Reactor fuel temperature 
Reactor exit lithium temperature 
Reactor inlet lithium temperature 
Argon temperature 


Time Final value 

of 

peak, °R 
sec 


Reactivity step, 1 cent 


11.2 6.2 103 

9.6 5.3 122 

8.0 4.4 153 

4.2 


Argon flow rate, 5 percent 



Reactor fuel temperature 

-2. 1 

-1.2 

34 

Reactor exit lithium temperature 

-4.0 

-2.2 

51 

Reactor inlet lithium temperature 

-7.6 

-4.2 

78 

Argon temperature 

-6.2 

-3.4 

61 



Lithium flow rate step, 25 percent 


Reactor fuel temperature 

-6.8 

-3.8 

3. 1 

Reactor exit lithium temperature 

-14. 1 

-7.8 

2.6 

Reactor inlet lithium temperature 

13.4 

7.4 

12.0 

Argon temperature 

-3.2 

-1.8 

2.0 



Transfer function approximation 


Peak over- Time Final value 
shoot of 

peak, 
sec 







2.7 

-3.2 

-1. 8 

29 

6.4 

3. 6 

.7 

-5.5 

-3. 1 

42 

2.0 

1. 1 

-1.9 

-9.0 

-5.0 

57 

-2.6 

-1. 4 

-1.2 

-7. 6 

-4.2 

42 

-1. 7 

-.9 


-.7 

-10.9 

-6. 1 

3.8 

-8.6 

-4.8 

-3.9 

-12.7 

-7. 1 

1.9 

-8.5 

-4. 7 

7.3 

11.9 

6.6 

8.0 

11. 5 

6.4 

-.6 

-7.5 

-4.2 

1.9 

-4. 1 

-2. 3 


Lithium flow 


Gas power 
conversion loop 


Heat 

exchanger 


1 U Argon flow 

Figure 1. - Primary loop of nuclear Brayton space powerplant. 




















































Lithium flow 


14.8 in. 
(37.6 cm) 



Void; 0.202 in. (0.513 cm) diam 
Fuel; 0.624 in. (1.58 cm) diam 
Ciad;0.75 in. (1.91 cm) diam 

Lithium cooiant passage; 

0.83 in. (2.11 cm) diam 

Segment 1 


- Segment 2 


■ Segment 3 


Figure 2. - Singiefuei pin model of reactor core (not to scale). 


Argon flow 
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Temperature change, 




0 100 200 300 400 500 0 100 200 300 400 500 

Time, sec 

(b) One delayed neutron group and three axial segment (d) One delayed neutron group and one axial segment with 

for reactor and heat exchanger. delta temperature of lithium to wall assumed to be zero. 



S S (e) Linear approximation of simplified equations. 


ev 



Time, sec 


(f) Inverse Laplace transform. 

Figure 4. - Time response about design operating point of reactor exit lithium temperature to 1-cent step disturbance in reactivity. 
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Temperature change, 





oH 




(a) Original simulation. 


(c) One delayed neutron group and one axial segment for 
reactor and heat exchanger. 




Time, sec 


(b) One delayed neutron group and three axial segments (d) One delayed neutron group and one axial segment with 

for reactor and heat exchanger. delta temperature of lithium to wall assumed to be zero. 


Figure 5. - Time response about design operating point of reactor exit lithium temperature to 5-percent step disturbance in argon flow rate. 






Temperature change, 





Figure 6. - Concluded, 
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